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Abstract: This is a second paper in our ongoing calculation of the next-to-next-to-leading 
order (NNLO) QCD correction to the total inclusive top-pair production cross-section at 
hadron colliders. In this paper we calculate the reaction qq — > ti+qq which was not considered 
in our previous work on qq — > ti + X [1] due to its phenomenologically negligible size. We 
also calculate all remaining fermion-pair-initiated partonic channels qq', qq' and qq that 
contribute to top-pair production starting from NNLO. The contributions of these reactions 
to the total cross-section for top-pair production at the Tevatron and LHC are small, at the 
permil level. The most interesting feature of these reactions is their characteristic logarithmic 
rise in the high energy limit. We compute the constant term in the leading power behavior 
in this limit, and achieve precision that is an order of magnitude better than the precision of 
a recent theoretical prediction for this constant. All four partonic reactions computed in this 
paper are included in our numerical program Top++. The calculation of the NNLO corrections 
to the two remaining partonic reactions, qg —> ti + X and gg — )■ ti + X , is ongoing. 
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1 Introduction 

Until very recently top-pair production at hadron colliders was analyzed in improved next-to- 
leading order (NLO) QCD. Broadly speaking, one can identify two such approaches, usually 
referred to in the literature as resummed and approximate NNLO, the latter simply being the 
truncation of the former to order 0{a%). 

The improved-NLO approximation to the NNLO cross-section is based on the next-to- 
next-to-leading log (NNLL) threshold approximation [2, 3] and also includes Coulombic terms 
[4] through NNLO. This approach is valid close to absolute threshold and the approximate 
results are added to the well known NLO [5-7] and NLL [8] results. Alternatively, in Refs. [9- 
11] the resummed NNLL (and, by truncation, the approximate NNLO) total inclusive cross- 
section was derived from the resummed differential one. 

A number of phenomenological studies have been presented in the literature [9-17]. Crit- 
ical comparisons of the various approaches can be found in Refs. [15, 17, 18]. 

As was demonstrated in Ref. [17], a true improvement in the theoretical precision in 
top-pair production at both Tevatron and LHC can be expected only upon inclusion of the 
full NNLO correction to the partonic cross-section. In our recent paper [1] we computed the 
dominant correction to top-pair production at the Tevatron confirming the expectations set 



At present theory agrees with data from the Tevatron and LHC [19-30] and the NNLO 
theoretical prediction for the Tevatron [1] has significantly smaller uncertainty than the 



in Ref. [17]. 
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existing experimental measurements. Theoretical predictions [31] based on the so-called 
BLM/PMC approach [32-34] have recently appeared and they, too, exhibit very small theo- 
retical uncertainty. A detailed comparison with the exact NNLO result is of interest and will 
be performed elsewhere. 

In Ref. [1] we argued that the purely fermionic channels that were not included there were 
phenomenologically negligible. Nevertheless, these reactions exhibit a logarithmic rise at high 
energy which could make them relevant for the description of lighter quarks, like charm or 
bottom, or for top-pair production at future higher energy hadron colliders. To that end, 
in this work we complete the calculation of all fermion-pair initiated contributions to top- 
pair production at hadron colliders through NNLO. Specifically, we calculate the reactions 
qq — » ti+ X, qq' — > tt + X and qq' — > tt + X, with q' ^ q. We also complete the computation 
of the qq initiated reaction [1] by deriving the result for the numerically subdominant reaction 
qq — )• tt + qq which was not considered in Ref. [1] . The results derived in this paper fully 
confirm our expectations [1] about the size of these reactions. Within the numerical accuracy 
of our calculation, we confirm the exactly predicted [35] leading high-energy logarithmic term 
of these reactions. Moreover, we are able to extract the subleading constant term in the high- 
energy expansion of the cross-sections. Our result is consistent with a very recent prediction 
[36] for this constant and improves the precision with which this constant is known by one 
order of magnitude. 

The paper is organized as follows: in section 2 we specify the reactions that we compute 
in the present work. Our calculational approach is explained in section 3. Factorization of 
collinear singularities is worked out in section 4, while in section 5 the scale dependence is 
derived. The explicit results are given in section 6. We conclude with section 7, where we 
discuss the phenomenological significance of the results computed in this work. 



2 Notation 

Following the notation of Ref. [1], the total inclusive top-pair production cross-section reads: 

^tot = V/ d/3$ ii (/3, Ai 2 ),T lj (/3,m 2 ,/i 2 ), (2.1) 
id Jo 

where i,j run over all possible initial sate partons, /3 max = a/1 — 4m 2 /S' with y/S the cm. 
energy of the hadron collider and f3 = \/l — 4m 2 / 's is the relative velocity of the final state 
top quarks having pole mass m and produced at partonic cm. energy yfs. 
The partonic flux appearing in Eq. (2.1) reads 

<M/^ 2 ) = 4* (t^T ■ ( 2 - 2 ) 

where, as usual, the partonic luminosity is 

dj(x,n 2 ) = x(fi ® fj) (x,fj?) = x I dy dz S(x - yz)f i (y)f j (z) . (2.3) 

Jo Jo 
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The scale \i in Eq. (2.1) stands for both the renormalization (hr) and factorization scales 
(Hf)- For = = H the NNLO partonic cross-section for the reaction ij — > ti + X reads 



+ 4 



^ ) + ^' 1) + ^ 2) 



(2.4) 



where L = In f/i 2 /m 2 ) and as is the MS coupling renormalized with iV^ = 5 active flavors at 
scale [i 2 . The functions erj^' depend only on /3. 

The partonic cross-sections are known exactly [5-7] through NLO. The scale controlling 

(2 1) (2 2) (1) 

functions and o\s can be easily computed from the NLO results <j\a , see section 5. 

The dependence on fj,R ^ \ip can be trivially restored in Eq. (2.4) by re-expressing ots(fJ*F) 
in powers of as(/J.R); see for example Ref. [12]. 

The dominant, phenomenologically relevant part of the reaction qq — )• ti + X was com- 
puted through NNLO in Ref. [1]. In this paper we compute the NNLO corrections to the 
remaining part of this reaction as well as the three new fermionic reactions qq —> tt + X, qq' — > 
tt + X, qq' ->tt + X, with q 1 / q: 



qq 



tt + qq\ m 
ti+qq , 
ti+qq' , 
ti+qq. 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



The label NS in Eq. (2.5) implies that we consider only those squared diagrams where 
the qq final state is not produced by a pure gluon splitting. Such (singlet) contributions are 
equal for all massless flavors, contribute with a power of Nl, and have been included in the 
calculation of Ref. [1]. 

The currently unknown contributions to ti production at NNLO are the qg and gg initi- 
ated reactions. They will be the subject of a future publication. 



3 The calculation 

All four partonic reactions (2.5,2.6,2.7,2.8) considered in this paper are of the so-called Double- 
Real type, i.e. they all have four-particle final states. They are computed with the STRIPPER 
approach of Refs. [37, 38]. We organize the calculation of the bare diagrams in the following 
way: we directly compute the bare contributions 

Ve{qq~' -> tt + qq') and a e {qq' ->• ti + qq') (3.1) 

to the reactions (2.6,2.7) while the two remaining partonic reactions (2.5,2.8) are computed 
as differences with respect to the reactions (2.6,2.7) 

af lS (qq) = a £ (qq -> tt + qq\ m ) - a £ (qq' -> ti+qq") , (3.2) 
af s (qq) = a £ (qq -> ti+qq) - a £ (qq' -> ti + qq') . (3.3) 
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The differences are derived from pure interference diagrams and vanish both at threshold 
(3 = and in the high-energy limit /3 = 1. 

The subscript e appearing in Eqs. (3.1,3.2,3.3) emphasizes that these are bare cross- 
sections, containing collinear singularities starting from 1/e 2 . To subtract these singularities 
and obtain the finite partonic cross-sections a, one needs to perform collinear factorization, 
which we describe next. 



4 Collinear factorization 

The description of the collinear factorization deserves some attention since for the reactions 
considered in this paper it has not been spelled out in the literature. Moreover, the collinear 
factorization for the reaction (2.5) represents a nonstandard contribution to the reaction 
qq — >■ tt + X and, for consistency, was suppressed in Ref. [1]. We take the opportunity to 
describe it in this work. 

In the notation of Eq. (2.4), and setting fj, = m, the bare partonic cross-sections read 

d ^ P ) = ^ + asvff (e, p) + o? s ~a^{e,p) + ...}. (4.1) 

They are defined in d = 4 — 2e dimensions and expressed in terms of the dimensionless variable 
p = Am 2 /s = 1 — (3 2 . To obtain the finite MS-subtr acted partonic cross-sections <Jij(p) one 
has to factor out the initial state collinear singularities: 1 



P 



E 

hi 



Oki{x) 



00- 



(4.2) 



defined as an expansion in (a,s/(27r)) n . Through NNLO we have: 



The MS collinear counterterms V are expressed through the space-like splitting functions 



Tij(e,x) = 5ij5(l - x) + a s T^\e,x) + a 2 s vf^(e,x) 

M £ >*) = -2^— r- ' 

2 



2e % 3 



(4.3) 



with [3 = IIC4/6 — Nl/3 and as the renormalized coupling at scale hr. 



1 We note a typo in Eq.(7) of Ref. [7], where a and a have been exchanged. This typo does not affect the 
rest of Ref. [7]. 
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It is more convenient to introduce the functions s) 4 



and sf^ denned as (s,p) 



a^\e, p)j p and sf^\p) = cr^j\p)/p. In terms of these functions the finite cross-sections read: 
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(4.4) 



,(i) 



For brevity, above we have suppressed the dependence on e and x. 

Next we consider the qq, qq' and qq' initiated reactions with q' ^ q. Introducing the 
notation q = (q,q',q'), we get: 
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99 ' 



(4.5) 



The function P , reads: 

99 



p(i) 



99 



P 



99 



+ P 



P 



(i),s 



'/'/ 



99 



if 
if 



(4.6) 



where P^' S , pW' y are the singlet and (gq) valence NLO splitting functions in the notation 
of Ref. [39]. The NLO functions sj^ appearing in this section are needed through order 
0(e 1 ). We have derived the subleading 0(e 1 ) terms by extending the results of Ref. [7]. For 
the manipulations involving harmonic polylogarithms [40] we have used our own software, as 
well as the program HPL [41]. All integral convolutions are computed numerically. 

Next we consider the qq initiated reaction. The NLO coefficient function s^j can be 

(2) 

found in Ref. [7]. The complete NNLO cross-section s„q reads: 



; ( 2 _) 
'99 
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b qq ' 
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99 
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99 
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e 99 



1 
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p(0) 
99 



p(°) + - f JL 

: V27T 



99 



g(l) 

b gg 



p(0) + 1 f J_ 



(4.7) 



where the NLO splitting function reads Pgg 



p(l),S p (i),v 
r gg r gg 



We recall that in Ref. [1] a subset of the Double-Real diagrams corresponding to the 
partonic process (2.5) were neglected due to their small size. Despite its phenomenologi- 
cal insignificance, however, the reaction (2.5) generates collinear singularities starting from 
1/e 2 which, in Ref. [1], were excluded from Eq. (4.7) in order to ensure consistent collinear 
subtraction. 
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The contribution to Eq. (4.7) that was excluded from Ref. [1] reads: 

2 



As (2) 

^ S <jr«?,NS 



+ - 



1 

1 



,(o) 
'<?<? 



J_p(0) (<> p(0) , I [p(l),S , c-(l) 
£ 2 9fl ^ 99 e L 39 99 



1 

2^ 



s(0) 

99 



p(0) 
^99 ' 



(4. 



i.e. the above result needs to be subtracted from the RHS of Eq. (4.7) to arrive at the result 
of Ref. [1]. 

The origin of the terms involving Pgq , Pqg^ and Pqq' S in Eq. (4.8) is easy to understand: 
they involve iterated emissions that are consistent with the initial and final states of the 



reaction (2.5). The only subtle contribution to Eq. (4.£ 



is the function Sqg that reads 



c(i) 



cl 



CfCa\ 
1 ) 



8-7x + 



2x 2 



In (a:) 



+ ■ 



1 + x' 



1 



7T 



y -21n(l 



x) ln(x) + ln 2 (x) — 2Li2(x) 



(4.9) 



The function Sqq^ is a partial contribution to the space-like splitting function Pqq' v and 
originates in the interference of the splitting process q — > q + q + q. This interference term 
cannot be extracted from Pqq^' V based on its color factor Cf{Cf — Ca/2) which is shared by 
a number of gluon emission diagrams that also contribute to Pqq V ■ We derive the function 
Sq]^ with the help of two independent direct calculations, which we describe next. 

First, by extending the results of Ref. [42], we compute directly the time-like (fragmen- 
tation) analogue Tqq of the function Sqq. Then, following Ref. [43], we analytically continue 
Tqq to space-like kinematics. For this particular contribution the analytical continuation is 
trivial and is just the usual replacement f(x) —> —xf(l/x) supplemented by standard an- 
alytical continuation across branch points for the involved logarithmic and polylogarithmic 
functions. Second, we identify the function Sqq as the second diagram from the class C in 
Fig. 7 of Ref. [44]. Since the result for this diagram is not available in that reference, we have 
directly computed it, following the methods of Ref. [44]. Both calculations lead to Eq. (4.9). 
Finally, as a by product of our calculation, we present for the first time the time-like 



(i),v 



function T, 



(i). 



11 



T« + S® = {CI-^ C 



ii 



ii 



) 



15(1 -x) + 7(1 + x) ln(x) + 1 + X ln 2 {x) 

1 — x 



(4.10) 



which has appeared in the literature on heavy flavor fragmentation [45, 46]. 

To derive the partonic reaction (2.5) we use the collinear subtraction term Eq.(4. 



s (2) 
?<2,NS 



s (2) + As {2) 



(4.11) 



Adding the result Eq. (4.11) derived in the present paper to the one derived in Ref. [1] we 
obtain the complete contribution to the qq initiated reaction qq — > ti+X at NNLO. Since the 
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contribution from the counterterm As^ NS (4.8) cancels in the complete qq — > ti+X result, the 
point-wise cancellation of the collinear singularities (within the numerical precision) observed 
both in this paper and in Ref. [1] serves as an additional check of our setup. 



5 Scale dependence 



The scale dependent terms a, 



' and a)^ 1 in Eq. (2.4) can be derived from: a) the 



rr (2 ' 1} and rr (2 ' 2) 

requirement that the measured hadronic cross-section o" tot in Eq. (2.1) be independent of 
the factorization scale [i through NNLO, b) the parton distribution functions fi satisfy the 
DGLAP evolution equations, and, c) the known running of the strong coupling constant. 



It is again natural to work in terms of the functions s 
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The powers of l/(2-7r) appearing in the above equations originate in the somewhat uncon- 
ventional choice of ag as the expansion parameter in Eq. (2.4). The expansion of the splitting 
functions is as in Eq. (4.3) where (3q is also defined. The two-loop beta-function coefficient 
reads ft = l7C\/6 - 5C A N L /6 - C F N L /2. 

The scale dependence for any specific reaction can be easily derived from the above 
equations. The expression for the qq reaction has been given in Ref. [12]. The scale-dependent 
terms for the reaction qq are not available in the literature and we give them here: 
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b qq 
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where the splitting function is given in Eq. (4.6) and s 



99 



99 
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'ij > 



(5.2) 



n 



0, 1 are the finite LO and 



NLO coefficient functions available in analytical form [7]. 

We have computed all convolutions numerically and produced our own fits for all scaling 
functions. We have implemented them in the program Top++ [47]: the ones for the complete 
qq reaction in version 1.2 and the ones for the qq, qq' and qq' reactions in version 1.3. 



6 Results 

We calculate the coefficient functions a\- for the reactions (2.5,2.6,2.7,2.8) numerically in a 
number of points on the interval /3 £ (0, 1). For short, we will sometimes refer to the set 
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of computed points and their numerical uncertainties as "data". Specifically, the qq' and 
qq' initiated reactions are computed in 80 points, with /3go = 0.999, as was also done in 
Ref. [1]. For the qq initiated reaction (2.5) we have added the point /? = 0.9999, while for the 
(/(/-initiated reaction (2.8) we have added two more points (3 = 0.99375 and /3 = 0.9999. 

The reason for including these additional points is to more accurately constrain the high- 
energy (3 — >■ 1 behavior of the numerically extracted partonic cross-sections. As is well known, 
the partonic reactions considered in this paper exhibit logarithmic rise at high energy due to 
diagrams where the top-pair is emitted in the f-channel. The leading behavior in the limit 
/3 — > 1 (or, equivalently, in the limit p — > 0) of the partonic cross-sections for all four reactions 
(2.5,2.6,2.7,2.8) is 

°ff^u hf \ p ^^Hp)+<* + 0{ P ). (6.1) 
The constant c\ has been predicted exactly in Ref. [35]. Its numerical value is 



Cl = -0.4768323995789214. (6.2) 

We have verified that for all four reactions (2.5,2.6,2.7,2.8) our numerical calculations 
(with unconstraint fits) return values for c\ that are within 2% from the exact result (6.2). 
Therefore, having verified the consistency of our calculation with the exactly predicted leading 
logarithmic term, in all subsequent fits we impose the exact value for the leading logarithmic 
term. This allows us to extract the constant cq with maximum precision, which turns out to be 
high enough to significantly improve the approximate prediction that has recently appeared 
in the literature [36], and to derive fits that are highly accurate even in the limit /3 — > 1. 

As we already anticipated in Ref. [1], and confirm in this paper, the contributions from 
the all-fermionic reactions to the total inclusive top-pair production cross-section at present 
hadron colliders are negligible. A more detailed analysis will be performed in section 7. 

6.1 qq' and qq' initiated reactions. 

The results for the partonic cross-sections for these two reactions read: 



qq' 



I QQ 



ci ln(p) - f3 2 exp (f qf ) , (6.3) 

-0.740572 - 31.2117/3 2 - 0.31495/3 3 + 15.8601/3 4 - 1.64639/3 5 + 18.9767/3 6 

+ ln 2 (p) (-3.16565p + 12.3828/9 2 ) + ln(p) (-19.6977p - 16.1386p 2 + 4.17707p 3 ) , 



qq 1 



1 QQ 



Cl ln(p)-/3 2 exp , (6.4) 
-0.740558 - 23.4518/3 2 - 0.193073/3 3 - 5.97215/3 4 - 0.541402/3 5 + 31.8227/3 6 
+ ln 2 (p) (-3.29162^ + 15.9932^) + ln(p) (-21.3725/3 - 11.1642p 2 + 8.64746p 3 ) . 



The constant c\ is given in Eq. (6.2). The analytical expressions in Eqs. (6.3,6.4) are derived 
as global fits of the set of 80 points we compute numerically. The data, and the corresponding 
fits, are plotted on Fig. 1. 



-8 - 



OS: 



0.006 
0.005 
0.004 
0.003 
0.002 
0.001 


-0.001 



ij=qqbar' 

ii-qq' 

























































































0.2 



0.4 



0.6 



0.8 



(2) 

Figure 1. The computed results, including numerical uncertainties, for the partonic cross-sections a qq , 

(2*) 

(blue) and u qq , (red). The discrete results, computed in 80 points, are overlaid with the corresponding 
analytical fits (see text). Both cross-sections diverge logarithmically in the limit (3 — > 1. 



As is evident from Fig. 1 these partonic cross-sections vanish at threshold (3 = and 
diverge logarithmically in the high-energy limit /3 — > 1. The quality of the fits (6.3,6.4) is 
very high for intermediate and large values of (3, i.e. in this region the precision of the results 
is restricted by the numerical precision of our numerical evaluation. Fitting in the region of 
small (3 turns out to be more problematic, however, since the two functions are as small as 
0(1O -10 ) over a sizable range of (3. In this range the distance between the fits (6.3,6.4) and 
the data is large compared to the size of the numerical uncertainty. However, the absolute 
size of the deviation data — fit is below 0(1O -6 ), i.e. the inaccuracy of the fits is completely 
immaterial for any foreseeable phenomenological application of these results. We find the 
simplicity of the analytical expressions in Eqs. (6.3,6.4) very appealing. Because of their very 
high absolute accuracy we have implemented them in the program Top++ [47]. 

From the fits we extract the following values for the constant cq: 

\ -2.5173 from <rJJ 



Co (from Eqs. (6.3,6.4)) = 
We note that the values of cq extracted from both reactions are compatible within the numer 



(6.5) 

2.5186 from a qq >, . 
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ical uncertainty. In the following we turn to the estimation of the uncertainty in the extracted 
value for cq. 

Due to the global nature of the fitting procedure, one might wonder how the quality of 
the fit at low /3 affects the quality of the fits in the phenomenologically more relevant region 
of large (3 and, in particular, the extraction of the divergent behavior at /3 — > 1. To estimate 
the robustness of the extracted high-energy behavior of the partonic cross-sections we derive 
second set of fits, with the aim of fitting both the small- and large-/3 regions. These alternative 
fits fit the data within the numerical uncertainties, except for the first few points, where they 
deviate from data within about 10 times the size of the numerical error. The absolute size of 
this deviation is 0(1O -10 ). The analytical form of these fits is much more cumbersome and 
the values of their parameters are highly tuned. For this reason we do not present this set 
of fits explicitly. Moreover, the phenomenological implication of the differences between the 
two sets of fits is completely immaterial. 

The only place where the difference between the two rather extreme fits plays a role is in 
the very large /3 behavior of the partonic cross-sections and in the extraction of the constant 
cq. Indeed, from this alternative set of fits we obtain 



Again, the extracted values of cq from the two reactions are compatible. We take the difference 
between the two types of fits, Eq. (6.5) and Eq. (6.6), as a measure of the uncertainty in the 
extraction of the constant cq from our calculation, which we estimate around 5%. 

Next we compare our result for the constant cq with the corresponding prediction of 
Ref. [36]. The value for Co predicted in Ref. [36] has substantial uncertainty, slightly above 
50%, and is predicted in the range (—1.4305,-2.43185). We see that our value for cq is 
consistent with the prediction of Ref. [36], albeit at the end of the uncertainty range quoted 
in that reference, and has an order of magnitude better precision. Overall, the agreement we 
find with the prediction of Ref. [36] (which was derived with completely different methods) 
is a non-trivial check for both setups. 

Finally, we would like to point out that the knowledge of the high-energy behavior (6.1) 
of the partonic cross-sections alone is insufficient for meaningful collider phenomenology. The 
reason for this is that the high-energy expansion of the partonic cross-sections is not well 
converging and thus not a good approximation outside the range of /3 ~ 1; it is only relevant 
for the description of heavy pair production at very large /3 which is not the case for top-pair 
production at the Tevatron and LHC. 

(2) 

To visualize this better, on Fig. 2 we plot the cross-section a J, and its high-energy 
leading-power approximation (6.1,6.5). It is easy to see that the behavior of the two functions 
is dramatically different, by three orders of magnitude or more, outside the narrow range 
j3 ~ 1. We note that Fig. 2 looks similarly for any one of the reactions (2.5,2.6,2.7,2.8). 
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Figure 2. Comparison of the partonic cross-sections cr^J (blue) and its leading power behavior in the 
high-energy limit Eq. (6. 1,6.5) (red). The blue curve, same as the blue cure on Fig. 1, appears to be 
zero on the scale of the red curve outside the narrow range /3 1. 



6.2 qq and qq(NS) initiated reactions. 

As emphasized in section 3, we compute the contributions to the reactions qq and qq(NS) as 
differences with respect to, respectively, the qq' and qq 1 processes. These differences, Acr^ NS 

(2) (2) (2) 

and Aa qq = a qq — a qq , , vanish in both the threshold and high-energy limits and read: 

+ A ^S!ns . (6-7) 
(1.53647/3 3 + 10.7411/3 4 ) p - 24.3298/3V + (-4.50719/3 3 + 15.4975/? 4 ) p 3 
+ (2.90068/3 3 - 4.98808/3 4 ) p 4 - 1.26644/3 20 ln(/3) 

+ ln 2 (p) (0.327143p - 10.7669p 2 ) + ln(p) (3.86236p - 21.332p 2 + 17.4705/? 3 ) , 

ag) = Cl ln(p)-^ 2 exp(/ g? ) , (6.8) 
f qq = -0.740558 - 22.8129/3 2 - 0.191648/3 3 - 6.58031/3 4 - 0.537669/3 5 + 31.7872/3 6 

+ ln 2 (p) (-3.25313p + 15.8988p 2 ) + ln(p) (-21.0783p - 10.8176p 2 + 8.64557p 3 ) . 



a 
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Figure 3. The computed results, including numerical uncertainties, for the interference partonic 



cross-sections A<7^ NS (blue) and Acr^ (red). The discrete results, computed in respectively 81 and 
82 points, are overlaid with the corresponding analytical fits (see text). Both results vanish in the 
limit (3 -> 1. 



The constant c\ is defined in Eq. (6.2). The data and the fits for the functions Acr^ Ng and 

(2) 

Acr qq are plotted on Fig. 3. 

(2) (2) 

The results for the functions AffJ- NS and Aa qq demonstrate that these functions vanish 

(2) 

in both limits p = 0, 1. Therefore, the high-energy behavior of the complete reactions cr qqNS 

(2) (2) (2) 

and (jqq is determined by the functions a qq/ and a qq , discussed in detail in section 6.1. 

The quality of the fits in Eqs. (6.7,6.8) is similar to the ones in section 6.1. The absolute 
quality of Eq. (6.7) is quite good, with absolute difference data — fit below O(10~ 7 ) for small 

(2) 

/3. The quality of the fit of Acr^- NS beyond the first seven lowest-/? points is dominated by 
the uncertainty in the numerical evaluation. 

For the gg-initiated reaction we have performed two fits. The fit in Eq. (6.8) is performed 

(2) 

for the total contribution a qq . Its absolute quality is also high, with absolute difference 
data — fit below O(10~ 8 ) for small /3. The quality of the fit of a q 2 q beyond the first ten 
lowest-/? points is dominated by the uncertainty in the numerical evaluation. From this fit 
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we extract the following value for the constant cq: 

c (from Eq. (6.8)) = -2.5196, (6.9) 
which is consistent with the ones extracted from Eqs. (6.3,6.4). 

(2) 

Finally we have performed a tighter, higher quality fit for the difference Aa qq . We do 
not present it explicitly here for the same reasons explained in section 6.1. 

7 Discussion 

In the present paper we calculate the NNLO corrections to total inclusive top-pair production 
at hadron colliders from the six-fermion partonic reactions (2.5,2.6,2.7,2.8). The results in 
this work, in particular, complete the calculation of the NNLO correction to the reaction 
qq — > tt + X [1] . The contributions from these reactions have been discussed in the recent 
literature [36, 48, 49]. 

As we already anticipated in Ref. [1], and confirm with our present calculation, the 
contributions from the all-fermionic reactions are phenomenologically insignificant for top-pair 
production at present hadron colliders like Tevatron and LHC. The numerical contribution of 
all four reactions (2.5,2.6,2.7,2.8) to the top-pair production cross-section at the Tevatron and 
LHC is presented in table 1. Specifically, we present separately the results for Acr^ws) an d 





Tevatron 


LHC 7 TeV 


LHC 8 TeV 


LHC 14 TeV 


Actons) [pb] 


-0.0020 


-0.0097 


-0.0124 


-0.0299 


°"<?<?,(NS) [P b ] 


-0.0009 


-0.0001 


0.0021 


0.0464 


Call [Pb] 


0.0003 


0.0970 


0.1504 


0.7885 


o"tot [pb] 


7.0056 


154.779 


220.761 


852.177 



Table 1. Contribution to the total top-pair inclusive cross-section due to the reactions computed 
in this work: due to the reaction (2.5) alone, Aa^^s) and a qq ,(nS)) an d due to all four reactions 
(2.5,2.6,2.7,2.8) combined, er a n. As a reference point, our pure fixed order prediction for <7 to t is also 
given. 

a qq,(NS) due to the reaction (2.5) as well as the combined effect cr a u due to all four reactions 
considered in this paper. As a point of reference we also present in table 1 the pure fixed 
order NNLO prediction cr to t for the total inclusive cross-section. 

The contributions from the reactions (2.5,2.6,2.7,2.8) are in the sub-permil range, both 
for central values and scale variation, for Tevatron and LHC at 7,8 and 14 TeV. The numbers 
in table 1 are computed in fixed order QCD with version 1.3 of the program Top++ [47] with 
default precision, rat = 173.3 GeV, central scales and MSTW2008nnlo68cl pdf set [50]. 

The results of the present paper might potentially be of interest for the description of 
lighter quark production (6 or c) or for top-pair production at possible future high-energy 
hadron colliders. Only in such cases, due to the partonic flux being peaked towards larger 
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values of /3, the high-energy rise of the reactions (2.5,2.6,2.7,2.8) might become phenomeno- 
logically relevant. 

We derive high-quality analytical fits for the partonic cross-sections in all four reactions 
(2.5,2.6,2.7,2.8). Our fits have the exact leading logarithmic behavior [35] in the high-energy 
limit. Therefore, any numerical difference in this limit due to the imprecision of our fits 
behaves no worse than a constant at large /3, i.e. as CQ Xact — Cg +0(p). Based on our findings 
in section 6.1 we estimate 

c exact _ c fit < O ( 10 -l) . (7.1) 

On the other hand, up to the point /3so = 0.999, our fits are quite accurate, typically much 
better than 1%, and thus a very good representation of the exact result. Therefore, it is only 
in the region beyond the point Pso = 0.999 where the difference (7.1) might start accumulating 
error. Barring extreme cases, however, we believe that all NNLO partonic cross-section fits 
derived by us so far are under good theoretical control in the full kinematical range. 

Finally, we would like to stress that the knowledge of the high-energy behavior (6.1) of 
the partonic cross-sections alone is insufficient for meaningful phenomenology. The reason for 
this is that the high-energy expansion of the partonic cross-sections is not well converging; see 
Fig. 2. The high-energy expansion is only relevant for the description of heavy pair production 
at very large /3 which is not the case of top-pair production at the Tevatron and LHC. 

Work on the calculation of the NNLO corrections to the two remaining partonic reactions 
qg tt + X and gg — > ti + X is ongoing and will be presented in a forthcoming publication. 
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